Transient entanglement in a spin chain stimulated by phase pulses. 
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Dynamics of the one-dimensional open Ising chain under influence of 7r-pulses is studied. It is 
shown that the application of a specific sequence of such instant kicks to selective spins stimulates 
arising of perfect dynamical pairwise entanglement between ends of the spin chain. Analytic formulas 
for the concurrence dynamics are derived. It is also shown that the time required to perfectly 
entangle the ends of the chains grows linearly with the number of spins in the chain. The final 
entangled state of the ending spins is always the same and does not depend on length the chain. 



I. INTRODUCTION 

Generation and distribution of entanglement between 
remote particles is a necessary requirement for quantum 
information transfer [l[ , Q . It turns out that spin chains 
are one of the most promising quantum channels for cre- 
ation of spin-spin entanglement Q. It appears, that 
in spite of existence of correlations in spin chains, the 
amount of entanglement generated between two particu- 
lar spins rapidly decays with the distance between them 
3, [J] even near the point of quantum phase transition 
5(. So that, the implementation of efficient mechanisms 
for entangling spatially separated subsystems is still a 
challenging problem, especially for massive particles. 

During the last years different types of spin-spin in- 
teractions have been extensively studied as possible can- 
didates to entangle remote particles. It has been proved 
that for certain types of spin-1/2 open Heisenberg chains, 
it is possible to entangle the ends of the chain by action 
of stationary Hamiltonians @ and implementing different 
mechanisms like mirror symmetry dimerized models 
or weak couplin gs a t the ends of the chain 0, perform- 
ing a quench [^|. [lpj] . or considering alternating couplings 
between the spins [ill ]. Also, several dynamical schemes, 
such as application of optimized time-dependent mag- 
netic field [12J and periodically time-dependent nearest- 
neighbor coupling [13| have been proposed to faithfully 
entangle the ends of an open spin chain. Unfortunately, 
only a partial (although a quite large) entanglement can 
be achieved with these schemes. Besides, owing math- 
ematical complexity of the problem, the main tool for 
studying such systems remains numerical simulations. 

The simplest type of interaction between massive par- 
ticles is the homogeneous Ising chain. Being the simplest 
type of interaction, it has been a candidate for imple- 
mentation of quantum information algorithms since the 



very beginning 14]. It worth noting a recent proposal 



for measurement-based quantum computation [la ] using 
homogeneous Ising-like interaction. Also, this kind of in- 
teraction accompanied by the global tt pulses allows a 
perfect qubit transport, quantum mirrors and universal 
quantum computation [161 ]. 

On the other hand, the possibility of application of 



pulses to quantum systems in order to preserve the co- 
herence was discussed several years ago [17| , and recently, 
the idea of applying optimized tt pulses at irregular inter- 
vals of time was proposed and studied by Uhrig [l8[ for 
efficient control of dephasing in spin systems and later 
applied for preserving entanglement in dephasing envi- 
ronments jl9j ]. 

In this article we propose a simple analytical scheme 
which allows to create a transient entanglement between 
the ends of an homogeneous spin chain by means of appli- 
cation of a sequence of tt pulses on some particular spins. 
We _apply tt pulses not to all qubits (as it was proposed 
in [16] to control a qubit transport) but only to a part of 
them, which results in a perfect transient entanglement 
between the ends the chain. 



II. THE MODEL 



Let us consider an open chain of N spins 1/2 with ho- 
mogeneous Ising-like interaction (the coupling constant 
is taken to be unity), governed by the Hamiltonian 



H 



N 



(1) 



and initially prepared in the coherent superposition 
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Since \^o) is not an eigenstate of the Hamiltonian 
Eq.(HJ), some specific spin-spin correlations arise during 
the Hamiltonian evolution, U(t) = e~ ltH . The evolution 
of the density matrix p can be found in the following 
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closed form, 

m = u(t) i* ) (* | uHt) 
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In the Hamiltonian ([1} only nearest neighbors interact, 
and a dynamical pairwise entanglement is generated be- 
tween them. As it can easily be seen from <j2j) > the entan- 
glement dynamics between all the pairs in the middle of 
the chain is the same, the concurrence [20| . is given by 
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max 
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where j = 2, . . . , N — 2 and its maximum is ~ 0.31. Since 
the spin chain is open, the concurrences for the spins 1 
and 2, and the spins N — 1 and N are different that given 
in Eq. ©, 



C12 — Cat-in — 



sini 
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reaching the maximum of 1/2 at t = (2n + 1)tt/2, with 
n = 0,1,.... 

It is important to stress that non neighboring pairs 
of spins will never get entangled under the Hamiltonian 
evolution. 



III. ENTANGLING THE ENDS OF THE CHAIN 
WITH tt-PULSES 

In order to generate entanglement between the ends 
of the chain let us apply instant 7r pulses to the first 
N — 1 spins at certain times tj. Such kicks correspond 
to rotations around the y-axis, and for the j-th spin are 
represented by the operator Rj — e~ l7TSvj / 2 . The main 
idea is as follows: initially all the spins are in eigenstates 
of s x j so that, in the course of the evolution governed 
by Eq.([T]), some spins become correlated. At the instant 
t\ = 7r we apply a 7r-rotation to the first N — 1 spins, 
producing flip of those spins. Afterwards, the evolution 
under the Hamiltonian Eq.(fT|) allowed to continue leading 
to the recoupling of all the spins. Then, at the instant 
t<i = 2tt we apply the inverse 7r-rotation, allowing the 
system to evolve further and repeat this cycle until N — 2 
spin flips are performed at times tj = jit, j = 1, 2, N — 
2, so that 

|*(t)) = U(t-Tr(N -2))R ( - 1)N ~ 1 U(n)... 



where R = ^Ufs^Rj. Using © we now proceed 
to compute the reduced density matrix and the (time- 
dependent) concurrence C(t) |20| . 

The three-spin case is special, since the middle spin is 
coupled to both ends. For a three-spin chain a single it 
pulse at t\ — it, applied to the first and second spins, 
produces entanglement between the ends. It is straight- 
forward to obtain the concurrence for the first and third 
spins: 



C 13 (t) = 



cos 




? £, t > 7T 
t < 7T 



i.e. the first and third spins, which were not entangled 
before the pulse (t < 7r), become completely entangled 
at t = 2nk, k = 1, 2, . . ., while C 12 {t) = C 23 (t) = for 

t > 7T. 

For longer chains the dynamical behavior of the pair- 
wise entanglement is rather unusual. After the first kick 
and before the (N — 2)-nd kick the pairwise entanglement 
between any pair of spins disappears. Since the system is 
closed, the pairwise entanglement transforms into a mul- 
tipartite entanglement during that time. Moreover, after 
the (N — 2) — nd kick it is possible to obtain the following 
analytic expression for the concurrence between the ends 
of the chain: 

a) for the even number of spins 



C(t) 



and 



max 



0. 



b) for the odd number of spins 



C(t) 



0, I COS : 



0, 



t> (N — 2)ir 
t<(N- 2)tt 



t>(N- 2)ir 
t<(N- 2)tt 



(5) 



(6) 

while the nearest neighbors in the middle of the chain 
recover their original transient entanglement Cjj+t(t) = 
max[0, | sint| - sin 2 t/2}/2 for times t> (N — 2)tt. 

Thus, the subsystem containing the first and the last 
spins of the chain stays in an unentangled state un- 
til the last kick is applied. As a result the subsys- 
tem reaches a maximally entangled state at the mo- 
ments t — (N + 2k— l)7r, k = 0,1,... , being unen- 
tangled with the rest of the chain, in accordance with 
the monogamy property [2l|. So that, the time needed 
to completely entangle the ends of the chain linearly de- 
pends on the length of the chain. In particular, at the 
instant t — (N — l)ir the state of the chain is as follows: 



I*) 



1 



2 (JV-l)/2 



JV-1 

n[i 

J=2 



[|00) 1Ar + |ll) 1Ar + i(|01) 1JV + |10) 1JV )]. 

On the other hand, the concurrence between the spins 
1 and 2, and the spins N —1 and A, remains zero for any 

t > 7T. 



3 



0.8 
0.6 
0.4 
0.2 


I 

1 

0.8 
0.6 
0.4 
0.2 




\ 










1 (a) 










I 








\ - 


i 










I 










i 

1 I 










i 1 




ii i/ 


\ '/ 


I 


i II 








li 


4 - '/ 


r 




V '/ 


It 























10 



15 20 25 30 



\ 1 1 








» (b) 








I 








i 
i 








. i 

1 
i 


7 


H '/ 


I 1 1 


\ 


'/ 


It '/ 


V 1 


\ 


if 


It if 


V i 




/ | 






' V t 2 t 3 : 
1 ' I i 


<4 \ / 








15 20 25 30 

t 



FIG. 2: (Color online) The concurrence C for chain of 4 spins 
as a function of the moments of applied kicks at fixed time 
t = 3tt. 



FIG. 1: (Color online) Dynamics of concurrence C(t) for the 
chain with (a) 4 spins for ti — n, ti — 2n; and (b) 7 spins 
for ti = 7r, t2 — 2n,ts = 37r, ti = 47r, t$ = 5-7T. Dashed line 
shows Tr{p\ N ), arrows indicate when the kicked pulsed are 
applpied. 

In Fig. (JT|) we show the dynamics of the concurrence 
C(t) for the chain with (a) 4 spins and 7r-pulses applied 
to the first 3 spins at tj — jir, j — 1, 2, and (b) 7 spins 
and 7r-pulses applied to the first 6 spins at tj = jn, 
j = 1, ...5. The dashed line shows the purity of the of 
subsystem composed by the first and the last spins of 
the corresponding chain, Tr(p\ N ). Observe that, due to 
the Hamiltonian evolution, the concurrence and purity 
keeps oscillating after application of the last kick. 

In order to study the behavior of concurrence C(t) 
when the times tj of application of the 7r-pulses are not 
commensurable, we have numerically calculated C(t) = 
C{t,t\,t2) as a function of t\ and ti > t\ at different 
fixed times. In FigJ5]we plot the concurrence C(t) for 4 
spins at time t = 37r (compare with Fig[T] (a)) when t± 
G [0.1,5] and ti £ [5.1,9] . We can observe that C{t\,t2) 
has a pronounced maximum, C = 1, at t± — n, and 
ti = 2ir and smoothly decreases when t\ t 2 deviate from 
these optimal values. 

IV. DISCUSSION 

We have studied the dynamics a one-dimensional open 
Ising chain of N spins assisted by phase kicks at times 



tj = jn, j = 1,2, ...,7V — 2. It is found that the ap- 
plication of N — 2 instant pulses to N — 1 spins leads 
to arising of transient perfect entanglement between the 
first and the last spins of the chain. One interesting re- 
sults is that, if the number of pulses is less than required, 
then every pairwise concurrence would be zero. The pe- 
riodic behavior of the concurrence between the ends of 
the chain, Cuv(t) reaches the maximum value C = 1 at 
some specific times even when the pulses are stopped be- 
ing applied and the chain evolves under the Hamiltonian 
(UJ) only. This effect substantially enlarges possible ap- 
plications of the low-dimensional spin chains in quantum 
information technology. 

On the other hand, in this particular model we may 
observe an important example of a non-trivial effect of 
local transformations on the entanglement production in 
non-linear systems. 

It is easy to see that combining the present scheme 
with quench one can entangle any two spins r and s, 
1 < r < s < N, in the chain, i.e. to model a quantum 
router. Really, we just have to "disconnect" the chain 
[Io| from 1 to r — 1 and from s + 1 to N and then apply 
the above discussed sequence of pulses. 

It worth noting that similar local transformations in 
the form of instant pulses were also used for another pur- 
poses: for atomic squeezing enhancement in the Dicke 
states [22^ . for intensification of the entanglement in 
continuous-variables [23| and in two-spin [24j systems. 

This work is partially supported by the Grant 106525 
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